Relaxation dynamics in fluids of platelike colloidal particles 
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The relaxation dynamics of a model fluid of platelike colloidal particles is investigated by means of 
a phenomenological dynamic density functional theory. The model fluid approximates the particles 
within the Zwanzig model of restricted orientations. The driving force for time-dependence is 
expressed completely by gradients of the local chemical potential which in turn is derived from a 
density functional — hydrodynamic interactions are not taken into account. These approximations 
are expected to lead to qualitatively reliable results for low densities as those within the isotropic- 
nematic two-phase region. The formalism is applied to model an initially spatially homogeneous 
stable or metastable isotropic fluid which is perturbed by switching a two-dimensional array of 
Gaussian laser beams. Switching on the laser beams leads to an accumulation of colloidal particles 
in the beam centers. If the initial chemical potential and the laser power are large enough a preferred 
orientation of particles occurs breaking the symmetry of the laser potential. After switching off 
the laser beams again the system can follow different relaxation paths: It either relaxes back to 
the homogeneous isotropic state or it forms an approximately elliptical high-density core which is 
elongated perpendicular to the dominating orientation in order to minimize the surface free energy. 
For large supersaturations of the initial isotropic fluid the high-density cores of neighboring laser 
beams of the two-dimensional array merge into complex superstructures. 

PACS numbers: 61.20.Lc, 64.60.My, 64.70.Md 
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I. INTRODUCTION 

Fluids of platelike colloidal particles, e.g., clay sus- 
pensions, are of enormous fundamental and technolog- 
ical relevance because of their abundance and their dis- 
tinct properties due to the orientational degrees of free- 
dom of the constituting particles. These systems ex- 
hibit many interesting phenomena such as flocculation, 
gelation^ging, and even li quid crystal phase transitions 
III, 0, i, EH i 0, 0, 0, During the last decade 

quite some progress has been made in both the synthesis 
flSj and the theoretical description of the equilibrium 
structure [11, [11 [11, [13, H [11, S [111, [11 [SI, M IIBI of 
suspensions of platelike colloidal particles. However, un- 
derstanding the non- equilibrium properties of these sys- 
tems is still a big challenge. 

Systems out of equilibrium are commonly analyzed ei- 
ther by investigating the relaxation of the system into 
equilibrium after an instant change of the conditions, e.g., 
due to switching an external field, or by studying the sys- 
tem in a (non-equilibrium) stationary state. The present 
theoretical work is devoted to the former case of relax- 
ation within a model fluid of platelike colloidal particles; 
stationary states will be addressed in the future. 

In view of the spatial inhomogeneities expected to be 
found, which in equilibrium systems are most adequ ately 
described by density functional theory (DFT) [26l [27l. [28| . 
the current investigation is performed within the frame- 
work of a phenomenological dynamic density functional 
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theory (DDFT), which proposes an equation of motion 
for the particle number density profiles [1^. The latter 
are assumed to describe the state of the system com- 
pletely, as within DFT any relevant quantity is a func- 
tional of the densities. Moreover, DFT is reproduced as 
the stationary limit of DDFT. 

The motivation of the DDFT equation within this work 
follows the traditional reasoning known from treatments 
of time-dependent Landau- Ginzburg and Cahn-Hillard 
models in studies of critical dynamics, spinodal decom- 
position, and crystal growth [H [HI, [3^. |33.]: the larger 
the thermodynamic distance from equilibrium, the faster 
the change of the state of the system. In the present case 
the state can change due to translation as well as due 
to rotation of the platelike colloidal particles. Moreover, 
the total number of particles in the system is conserved 
whereas the orientation is not. Hence the DDFT to be 
detailed in Sec. [IT] is analogous to "model C" in the clas- 
sification of Hohenberg and Halperin [s^l ■ Alternatively, 
the present DDFT can be regarded as an elaborate phase 
field theory [35!, |36j with the order parameter tensor as 
the natural phase field. 

In recent years DDFT equations have been derived 
based on (overdamped) Langevin dynamics by approx- 
imating the time-dependent two-particle density [stI . [ssI . 
[bqI . [40| . As Langevin dynamics is considered a reason- 
able description for dilute colloidal suspensions, and as 
the isotropic-anisotropic liquid-crystal phase transitions 
in fluids of highly anisotropic colloidal particles take place 
at small number densities, DDFT is expected to be ap- 
plicable within the isotropic-nematic two-phase region of 
fluids of platelike colloidal particles. 

Here relaxation dynamics is investigated by consider- 
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ing a two-step switching process: First an initially homo- 
geneous stable or metastable isotropic state is brought 
out of equilibrium by switching on a two-dimensional ar- 
ray of Gaussian laser beams. The laser potential acts 
as an external potential which tends to form an inho- 
mogeneous equilibrium state. In Sec. IIIII the relaxation 
towards this inhomogeneous equilibrium state is traced 
by integrating the DDFT equation for the presence of 
the laser potential. After equilibrating the system in the 
presence of a laser potential, the laser beams are switched 
off and the dynamics is described by the DDFT equa- 
tion for vanishing external field. The system either re- 
laxes back to the initial isotropic state or it evolves into 
an anisotropic state. The various relaxation paths will be 
analyzed in Sec. lIVI This switching process has been cho- 
sen because it is expected to be realizable in experiments 
and, at the same time, numerically convenient boundary 
conditions can be used. As an illustration and in order to 
allow for quantitative comparison the model parameters 
have been chosen to describe an aqueous suspension of 
sterically stabilized gibbsite platelets (ilj . 

SectionlVl concludes the present work with a discussion 
of the applied formalism and the numerical results. 



II. FORMALISM 

A. Model fluid 

The system under consideration is a colloidal suspen- 
sion of monodisperse hard platclike particles dispersed 
within a continuous solvent. The colloidal particles are 
modeled by square cuboids which can take only one out of 
three mutual ly p erpendicular orientations (Zwanzig ap- 
proximation chosen to be parallel to the Cartesian 
axes. A particle is called a i-particle if its square face nor- 
mal points along the i-axis, i £ {x, y, z}. The local num- 
ber density of i-particles at position r is denoted by Pi(r) 
and the abbreviation g := (g^, gy, g^) is used for conve- 
nience. The structure of the model fluid is adequately 
described in terms of the total density g := gi and 
the order parameter tensor Q [43l | which, within Zwanzig 
models, is given by the diagonal form 



Here Sui denotes the Kronecker-(5. 

B. Equilibrium density functional theory 



(1) 



Density functional theory (DFT) [26!,l27| is the method 
of choice in order to investigate equilibrium properties 
of spatially inhomogeneous fluids [11]. An accurate 
version of DFT for the Zwanzig model, which repro- 
duces the exact second and third virial coefficients and 
which possesses the property of dimensional crossover, 



is the fundamental measure theory due to Cuesta and 
Martinez- Raton 4^ 45 L It has been applied to investi- 
gate monodisperse |l8l.l46j and polydisperse [i^l Zwanzig 
models. Within this framework, the free energy func- 
tional is given by 



me] 



J d^r (^^ft(r) (ln(ft(r)A3) - 1 + mr)) 



$(n(r)) , 



(2) 



where (3 is the inverse temperature, A denotes the ther- 
mal de Broglie wavelength, Vi (r) represents the external 
potential exerted onto i-particles at position r, and 



$(n(r)) = no(r)ln(l-ri3(r)) 



1 



n"2g(r) 



l-n3(r) +(l-n3(r))2 

describes the local excess free energy density. The latter 
is a function of the weighted densities na{r) := J2i ^a,i ^ 



gi{r) with i 
defined by 



denoting convolution and weight functions 



wiz,i(r 



0-{fxj Sxi)a.{ry , Syi)a[rz , Szi), 
b{rx: Sxi)a{ry , Syi)a,{rz , Szi), 
a{rx, Sxi)b{ry , Syi)a(rz, Szi), 
a{rx, Sxi)a{ry , Syi)h{rz, Szi), 
a{rx,Sxi)b{ry, Syi)b{rz, Szi), 
b{rx,Sxi)a{ry, Syi)b{rz, Szi), 
bijx, Sxi)b{ry, Syi)a(rz, Szi), 
bijx, Sxi)b{ry, Syi)b(rz, Szi), 



(4) 



where the abbreviations a{r,S) :— ^{S{§ + r) + 



r)) and b{r, S) 
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are used and S, 



qi 



s _ 

2 

denotes 



the extension of i-particles along the q-axis. The term 
" fundamental measure theory" is related to the geometric 
interpretation of spatial integrals of the weight functions 
Eq. (Ill) as particle number, linear extension in x-, y-, and 
z-direction, cross-sectional area perpendicular to the x-, 
y-, and z-axis, as well as particle volume, respectively. 

Within a canonical system the equilibrium density pro- 
files g'^'^ minimize F[g\ under the constraint 



J ~ const. 



(5) 



The corresponding Lagrange multiplier is the chemical 
potential fj,. With the local chemical potential 

l^^M-.^j^^, (6) 

g'"^ fulfills the Euler- Lagrange equation 

t,,{r,[g^^]) ^ ^l, (7) 

i.e., equilibrium density profiles render the local chemical 
potential as a function of position (r) and orientation (i) 
into a constant. 
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C. Dynamic density functional theory 

A system initially out of equilibrium will be driven to- 
wards equilibrium. Motivated by equilibrium DFT it will 
be assumed that the state of the system is described by 
time- dependent density profiles g{-,t). This assumption 
implies the neglect of hydrodynamic interactions, which 
depend on the velocity field. It is found that hydrody- 
namic interactions become more and more relevant for 
increasing packing fractions [H, [4^. As the proposed 
theory is intended to be applied to dilute colloidal sus- 
pensions, taking hydrodynamic interactions not into ac- 
count is considered to be a reasonable approximation. 

In order to investigate the temporal evolution of the 
system an equation of motion for g{-,t) has to be pre- 
scribed. Here the following form is proposed: 



dgi{r,t) _ I dg,{r,t) 
dt ^ \ dt 



dgtjr, t) 
dt 



(8) 



where the first term on the right-hand side describes the 
contribution due to the translation of particles with fixed 
orientation and the second term represents the contribu- 
tion due to the rotation of particles keeping the local 
total density g{r,t) constant. 

The translational part fulfills the continuity equation 



f dg,{r,t) 
\ dt 



E 



(9) 



where jiq describes the translational current of i-particles 
in ^-direction. Following Dicterich et al. f29| , the current 
is assumed to be proportional to the local density and to 
the negative gradient of the local chemical potential: 

3^g{r,[g{-,t)]) := -rigg,{r,t) . (10) 

OTq 

Here the translational diffusion matrix T with 



Til ,i = q 



(11) 



accounts for different diffusivity of platelike colloidal par- 
ticles parallel (Fn) and perpendicular (F^) to the sym- 
metry axis. 

The rotational part of Eq. ([8]) is modeled by a master 
equation 



dgi{r,t) \ ^ / dgi{r,t) 
dt Aot^V dt 

with the "gain" term 



^g»(r,^) 
dt 



loss 



(12) 



dgi{r,t) 
dt 



gam 



(13) 



J2 {f3^^^' (r, [g{-,t)]) - (3fi, {r, [g{- , t)])^ g,, {r , t) 



and the "loss" term 

loss 



( dgi{Y,t) 



\ dt 

- (r, [g{-,t)]) - (r, [^(^ t)])) g»(r, t) 



(14) 



Hence 

f dgdjc^y 
V dt 



i' 

(/3/i,(r,[^(-,<)]) -/3/i,.(r,[^(.,i)])). 



The proportionality factors of Eqs. p3)l and ([14]) equal 
the rotational diffusion coefhcient, which in terms of the 
rotational relaxation time r is given by 

In the limit of infinitely thin platelikc colloidal particles 
Brenner [50] has worked out translational and rotational 
diffusion coefficients. These expressions in the present 
notation lead to the rotational relaxation time 



(16) 



where 77 is the viscosity of the solvent and D denotes the 
diameter of the platelike colloidal particles, and to 



" 36 ' 



r ^ 



(17) 



Equations (gl), dH), (HU]), and ^ together with def- 
inition Eq. ([6]) and the density functional Eqs. ©-([I]) 
completely specify the dynamic density functional theory 
(DDFT) for Zwanzig particles to be investigated within 
this work. 

Integration of the equation of motion Eq. ([8|) with an 
initial configuration g{-,0) leads to the time-dependent 
density profile g which contains all spatial and temporal 
information. 

The proposed DDFT is consistent with equilibrium 
DFT as any equilibrium state g'^'^ fulfills Eq. ([7|) and 
therefore does not change under the dynamics repre- 
sented by Eqs. P. ((9)). (fTO)) . and p5)) . 



D. External potential 

Section IIIII studies the relaxation of the model fiuid 
within an external field. In view of conceivable exper- 
imental realizations this potential is chosen to model a 
two-dimensional array of Gaussian laser beams. Follow- 
ing the notation of Ref. 51] the force density onto an 
inhomogencous dielectric fluid is given by 



f(r) 



1 



-eoE2(r)Ve(r), 



(18) 



where Sq is the permeability of the vacuum, denotes 
the temporally averaged square of the local electric field, 
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and e(r) is the local relative permeability. A single Gaus- 
sian beam of width w and power P propagating parallel 
to the z-axis is described by [5l| 

E2(r) = 2 exp ^ (19) 

with the refractive index n and the speed of light c. The 
total force Fi(r) exerted on a single i-particle at posi- 
tion r is given by an integration of f over the particle 
volume. Given the refractive indices of the solvent risoiv 
and of the colloidal particle tIcoU, £ can be assumed to 
discontinuously change from nl^y^ outside the particle to 
"coll inside the particle. The aforementioned volume inte- 
gration hence reduces to a surface integration which can 
readily be performed in the case of Zwanzig particles. It 
is found that Fi(r) = — VVi(r) with the laser potential 
of a single beam 

V.{v) = - ^^"^°"~J'^°'^^^-'' ^(r.,^..)^(^y,g.O, (20) 

where Sqi denotes again the extension of i-particles along 
the g-axis and 

„(.5):.erf(^(,.+ f))-„t(:^(.-f))^,21) 

For ricou > n-soiv the laser potential is attractive. 

E. Choice of parameters and numerical method 

In order to obtain solutions of the DDFT equation 
specified in Sec. Ill CI one has to fix the model param- 
eters describing the particle geometry, the particle dif- 
fusivity, and the laser potential. In the following these 
parameters have been arbitrarily chosen to model the 
aqueous suspension (viscosity "q = 8.9 • 10~^Pa • s, refrac- 
tive index nsoiv — 1.33) of sterically stabilized gibbsite 
platelets (refractive index ricou = 1.58) AlOP of Ref. (4lj | 
(diameter D = 165nm, aspect ratio 12) irradiated by 
a two-dimensional square array of laser beams each of 
power P ~ lOmW and width w — 5D propagating in 
z-direction. Neighboring laser beams within the two- 
dimensional array are chosen to be 301? apart. 

Due to the periodicity of the total laser potential the 
actually infinite system is described in terms of a finite 
system with periodic boundary conditions under the in- 
fluence of one single laser beam centered at the origin. 
The density profiles, which vary only in x- and y-direction 
but not in z-direction, are defined on a square of side 
length 301? with periodic boundaries. Due to symmetries 
of the external potential only one quarter has actually to 
be coded which is done by means of a grid with spac- 
ing -J. This grid leads to slightly broadened isotropic- 
nematic interfaces as compared to calculations with finer 
grids. However, such a minute inaccuracy is considered 
irrelevant as compared to the numerical advantage gained 
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FIG. 1: a) The orientationally averaged laser potential V{x, y) 
exhibits an approximately axial symmetric attraction of plate- 
like colloidal particles by the laser beam, b) Deviations from 
the rotational symmetry appear on a smaller energy scale. 



by the relatively coarse grid. A particularly efficient nu- 
merical implementation is possible if one approximates 
the colloidal particles by infinitely thin platelets; an ex- 
ception is Eq. (PO)) where a thickness of ^ is used in order 
to avoid a vanishing laser potential. 

Figure [T] displays the orientationally averaged laser po- 
tential V :— ^ Vi as well as the local difference Vx — Vy 
of the potentials exerted onto x- and y-particles. The at- 
tractive laser potential is approximately axial symmetric 
on the energy scale f3~^ (Fig. [1^) whereas there are de- 
viations on the energy scale 10^^/3^^ (Fig- lb)- 

According to Eq. (|16p the rotational relaxation time 
is given by r « 214/is. The calculated translational dif- 
fusion coefficients (see Eq. PT]) ) compare well with the 
measured ones of Ref. [41[. The DDFT equation is in- 
tegrated with respect to time by means of the Euler- 
forward method with integration time steps of 

For the given set of parameters, one finds a first- 
order isotropic-nematic bulk phase transition at a re- 
duced chemical potential = -1.37271680846513 with 
/i* := (3fx — 31og(^). At this binodal (b) the isotropic 
bulk phase of density g'^°D^ « 1.144 coexists with the 
nematic bulk phase of density g'^°'^D^ « 1.576 and scalar 
order parameter 5 = (| cos(??)^ — i) « 0.83, where (•) 
denotes the thermal average and 19 is the an gle between 
the particle normal and the nematic director [43] . For il- 
lustrative purpose the reduced chemical potentials /i* are 

also expressed in terms of the supers aturation a :— ^ ^ia 
of an isotropic fiuid of reduced chemical potential /i*. 
The isotropic-nematic interfacial tensions for the nematic 
director pointing parallel and perpendicular to the inter- 
face normal are given by I3^\\D'^ — 2.051633048-10^'* and 
P^^D'^ = 3.832557464-10"'', respectively. The isotropic 
spinodal (is) is located at the reduced chemical poten- 
tial « —1.128 where the metastable isotropic state 
takes a density of Qf°D^ k, 1.244. The nematic spin- 
odal (ns) is located at the reduced chemical potential 
w —1.399 where the stable isotropic state takes a 
density of g'^°D^ « 1.133. 
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III. RELAXATION IN EXTERNAL FIELD 

This section describes the influence of the external 
laser potential Vj(r) (see Sec. Ill Dl) onto an initially ho- 
mogeneous isotropic fluid of platelike colloidal particles 
(see Sec. Ill A| ). Depending on the initial chemical poten- 
tial, two qualitatively different scenarios are possible. 

For reduced chemical potentials fi* < fi* the fluid stays 
isotropic throughout and the laser potential merely leads 
to an axial symmetric accumulation of colloidal parti- 
cles near the origin. For the parameters of Sec. Ill El the 
limiting reduced chemical potential is ji* ~ —1.6 (super- 
saturation a K, —0.080). 

However, if the initial reduced chemical potential ful- 
fills jjL* > pL* the initial isotropic symmetry of the fluid is 
broken. Initially homogeneous states within the range 
H* £ [fx* , exhibit the same qualitative behavior. 
Hence it is sufficient to describe the case of one initial 
state corresponding to, e.g., ^x* = — 1.25 (cr = 0.044). A 
representation of the solutions of the DDFT equation in 
terms of the total density proflle g and the order param- 
eter tensor component proflles Qxx and Qyy is displayed 
in Fig. H 

After a time oit = 50t slight distortions of the initial 
homogeneity due to the attraction of the colloidal parti- 
cles by the laser beam become visible in the total density 
proflle Q-, the density at the origin is p(0, 0)1)^ « 1.382. 
Orientational ordering, however, is still negligible. 

At time t = 290t more particles have been concen- 
trated near the origin {q{{),Q)D^ k, 1.450) at the ex- 
pense of the direct surrounding. Moreover, four max- 
ima of the total density profile g formed near the origin 
{Q{±l.bD,Q)D^ = (?(0,±1.5L»)£'3 ^ 2.250). The cubic 
symmetry of the total density profile g and the rectan- 
gular symmetry of the order parameter components Qxx 
and Qyy reflect the orientation dependence of the exter- 
nal laser potential Eq. ([^0]) (see Fig. [TJd). 

At < = 440r, the cubic symmetry of the total density 
proflle is broken by forming a bridge between two op- 
posite peaks present at time t = 290r. The density at 
the origin has increased to g(0, 0)D^ « 2.723. In the 
current case this bridge is formed by x-particles but a 
bridge made of y-particles is possible as well. It has 
been conflrmed that minute rounding-off errors in the 
very flrst integration step gives rise to the formation of 
bridges of one or the other orientation. The two peaks 
of y-particles, which do not form a bridge, have moved 
away from the origin and the local maximum of the to- 
tal density is given by £((±3.513, 0)D^ « 1.648. They are 
subsequently decreased because the nearby bridge forces 
the particles to take a;-orientation which afterwards join 
the bridge. 

At time t = 500t the two peaks of y-particles have dis- 
appeared completely. The approximately elliptical core 
(semiaxes of half height contour approximately 2.1D and 
2.6-D) of increased particle concentration near the center 
(p(0, Q)D^ K, 2.805) comprises mostly x-particles. In the 
following more and more particles from the surrounding 



are added to the high-density core giving rise to a de- 
crease of the total density outside the laser beam. 

From time t = 2000t onwards the system is almost 
equilibrated and the state practically does no longer 
change in time. The flnal high-density core has, as the 
orientationally averaged laser potential V (Fig. [2^), an al- 
most circular cross-section (radius of half height contour 
approximately 2.913). The total density at the origin 
is q{Q, 0)D^ « 3.530 whereas the minimum of approxi- 
mately 1.101 is found at a distance of about 8.513 from 
the origin. Due to the small flnite size (diameter 30Z3) 
of the considered part of the total system, the total den- 
sity at the (periodic) boundaries has been decreased to a 
value of g\oD^ « 1.104. 

The trajectory of the system for /i* > p,* could be sum- 
marized as follows: The laser potential, which is not axial 
symmetric (see Fig.[T]), attracts particles to the origin and 
thereby creates four nematic nuclei (two of x-particles 
and two of y-particles) which compete at the origin. If 
one pair of nuclei breaks the symmetry at the origin it 
forms a bridge whereas the other pair of nuclei decays. 
This scenario hinges on the creation of the nematic nu- 
clei and hence on a sufficiently strong laser potential: The 
limiting reduced chemical potential jl* is expected to de- 
crease upon increasing the laser power P. In particular, 
for sufficiently strong laser potentials, jl* can be located 
well within the isotropic phase of the bulk phase diagram, 
i.e., the symmetry breaking scenario can occur for initial 
states that correspond to a stable isotropic fluid. 

For a given initial reduced chemical potential fi* one 
can deflne a limiting laser power P{fJ.*) such that the 
symmetry breaking scenario will take place if and only if 
P > P{ii*). For the case n* = -1.25 (ct = 0.044) and 
the parameters of Sec. Ill El the limiting laser power has 
been determined approximately as P(— 1.25) « 3.5mW. 



IV. RELAXATION OF THE FREE FLUID 

The previous section described the equilibration path 
of an initially homogeneous isotropic fluid within a laser 
potential. The equilibrium structure in the presence of 
the laser potential has been attained after about t — 
2000t (see Fig. [5]). The present section is concerned with 
the temporal evolution of the fluid after switching off the 
laser beam at time t* :— 2000t. Three possible scenarios 
for the evolution of the free fluid have been identifled 
which are determined by the total density of the initial 
homogeneous isotropic fluid at time t = 0. 

If the total density at time t = is small {fi* ^ 
— 1.32, (T ^ 0.019), the high-density core present a.tt — t* 
(see Sec. lIII[) dissolves completely for i ^ cxd by releasing 
the excess amount of particles to the surrounding. Un- 
der these conditions the system relaxes back to the ini- 
tial homogeneous isotropic state. In the case of a super- 
saturated initial isotropic system {fi* € (/ij^, — 1.32), cr e 
(0, 0.019), see Sec. Ill Ep . classical nucleation theory would 
consider the perturbation due to the external potential 
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FIG. 2: Temporal evolution of an initially homogeneous supersaturated isotropic fluid of platelike colloidal particles under the 
influence of a laser potential in terms of the total density profile g and the order parameter tensor component profiles Qxx and 
Qyy. The length scale D and the time scale r correspond to the diameter and to the rotational relaxation time of the colloidal 
particles, respectively. The equilibrium densities at isotropic-nematic two-phase coexistence are given by gb°D^ ~ 1-144 for 
the isotropic state and by £iJJ°™_D^ ~ 1.576 for the nematic state. 



too small to cross the high free energy barrier between 
the weakly supersaturated metastable isotropic state and 
a stable nematic state. 

For intermediate initial supersaturation the high- 



density core reshapes with time attaining for t — > 
oo a finitely elongated, approximately elliptical shape. 
Figure [3] displays the structures of the almost equi- 
librated fluids for initial reduced chemical potentials 



7 




of constant area of the ellipse leads to the aspect ratio 
1.84249. This value agrees well with the value 1.85 of 
the equilibrium structures of Fig. [3] The calculation also 
shows that the aspect ratio depends only on the ratio of 
the interfacial tensions but not on the area of the ellipse. 



Within the terminology of classical nucleation theory 
the perturbation due to the laser potential has been 
large enough to cross the free energy barrier between 
the metastable isotropic state and a state containing an 
approximately elliptical nematic core surrounded by an 
isotropic fluid. In compliance with classical nucleation 
theory, the shape of the core is determined by a minimum 
of the surface free energy and the densities within and 
outside the core are approximately given by the isotropic- 
nematic coexistence values. However, whereas classical 
nucleation theory and its generalizations [s^l seeks for the 
fluctuation induced critical nucleus, which corresponds 
to the free energy " saddle point" between the metastable 
and the stable state, the perturbations considered in the 
present work are of external origin and the free energy 
barrier is crossed along "higher" paths. 



FIG. 3: Structure of a free fluid of platelike colloidal particles 
for reduced chemical potentials ^* £ { — 1.31, —1.27, —1.21} 
(supersaturations a € {0.022, 0.037, 0.058) at time t ^ t* + 
lOOOOr. The fluid has been exerted for a time interval t* = 
2000r to a laser potential (see Fig. [2| which afterwards has 
been switched off. These structures represent the equilibrium 
structures as they do not change in time anymore. 



/i* e { — 1.31, —1.27, —1.21} (supersaturations a E 
{0.022,0.037,0.058) at time t = t* + IOOOOt. These 
structures exhibit the same total density at the origin 
q{0, 0)D^ « 1.62, the same total density far from the ori- 
gin g\dD^ ~ 1.15, as well as the same aspect ratio 1.85 
of the half height contour. Note that the core extension 
of the equilibrium state can be as small as a few particle 
diameters (see case fi* = —1.31 {a = 0.022) in Fig. [3]). 
Moreover the high-density core imposes orientational or- 
der to such a degree that a preferred alignment of the 
platelike particles is found even within the dilution zone 
surrounding it. 

An estimate of the aspect ratio of the equilibrium half 
height contour can be obtained within a sharp-kink ap- 
proximation which describes the equilibrium structures of 
Fig. [3] by a nematic ellipse surrounded by isotropic fluid 
where the isotropic-nematic interface is described by a 
step function. Assuming an angle-dependent interfacial 
tension 7(1?) := 7|| cos('!^)^ -I- 7_l sin(i?)^, where d denotes 
the angle between the contour normal and the nematic 
director, and where 7|| and 7_l are the isotropic-nematic 
interfacial tensions given in Sec. Ill Ei one obtains an ap- 
proximate surface free energy i^surf by integration of 7 
along the contour of an ellipse. Minimizing Fgurf with re- 
spect to the semiaxes of the ellipse under the constraint 



If the initial reduced chemical potential fi* is suffi- 
ciently large the high-density core contains enough par- 
ticles to grow to the (periodic) boundary of the modeled 
part of the system where it merges with its " images" into 
a system-spanning superstructure. Figure [3] displays this 
behavior for the case /i* — —1.13 (cr = 0.087). Directly 
after switching off the laser beams at time t = t* a, circu- 
lar high-density core with total density q{Q, 0)D^ w 3.734 
is present. At time t ~ t* + lOOOOr the total density at 
the origin has decayed to g{0,0)D^ w 1.62 and the as- 
pect ratio of the half height contour is 1.81. These values 
agree with those of the cases displayed in Fig. [31 How- 
ever, the region of orientational order, which extends be- 
yond the dense core, has reached the periodic boundaries 
of the system. The high-density core thereby interacts 
with its "images" which leads to a further stretching as 
displayed for time t — t* -\- 15000t. Finally the high- 
density core merges with its "images" forming a ban- 
dlike structure with a total densities g{0,0)D^ « 1.58 
and g{ztl5D,0)D^ « 1.14. This structure corresponds 
to a lower free energy as compared to the approximately 
elliptical structure because the isotropic-nematic inter- 
face is now everywhere parallel to the nematic director. 
The time scale on which the high-density core stretches 
is much larger than the time scales of the processes ob- 
served in the presence of the laser potential (see Sec. lIIip 
because the chemical potential gradients and hence the 
currents are much smaller here. In an experimental real- 
ization the high-density cores at the centers of the indi- 
vidual laser beams will be oriented in different directions. 
Upon merging the high-density cores a complex network- 
like structure is expected to form. 



8 




Li? ^3 



I: 



i= i' + IOOOOt; 



-la -10 J a 10 la 

ilD 



-la -10 -a a lo la 




-13 -10 -3 3 10 15 



-IS -10 -a 3 10 IS 




-la -10-3 a 10 la 

i;'-D 




-la -10 -a a lo la 



FIG. 4: Growth of the high-density core within a free fluid of 
platelike colloidal particles for initial reduced chemical poten- 
tial ^* = —1.13 (supersaturation a = 0.087). The fluid has 
been exerted for a time interval t* = 2000r to a laser poten- 
tial (see Fig. [2]) which afterwards has been switched off. The 
growing high-density core ultimately merges with its "images" 
due to periodic boundary conditions. 



V. DISCUSSION AND SUMMARY 

The previous three sections introduced and applied a 
novel dynamic density functional theory (DDFT) for a 
model fluid of platelike colloidal particles. The particles 
are assumed to be infinitely thin and the orientational 
degrees of freedom are taken approximately into account 
by restricting the particle normal to three mutually per- 
pendicular directions (Zwanzig model [4^). Although 
qualitative, this model exhibits, in agreement with ex- 
perimental studies [Tol . [Tl| , an isotropic- nematic bulk 
phase transition at low number densities. Within the 
framework of DDFT the state is described by a set of 
density profiles and the proposed equation of motion is 
constructed in order to fulfill the conservation of the to- 
tal number of particles. Time dependence is generated 
by gradients of the local chemical potential which is de- 



rived from a density functional. The phenomenological 
parameters appearing in the equation of motion are de- 
termined by matching the dilute limit with single parti- 
cle translational and rotational diffusion 50]. This ap- 
proach does not take hydrodynamic interactions between 
colloidal particles into account because they are assumed 
to be small for the small densities within the isotropic- 
nematic two-phase coexistence region considered here. 

The proposed DDFT is applied to the investigation of 
relaxation dynamics of the model fluid. In a first step 
the relaxation of an initially homogeneous isotropic fluid 
in an external potential due to a two-dimensional array 
of laser beams (Fig. [1]) is studied. The restriction on the 
orientational degrees of freedom of the particles leads to 
structures of rectangular symmetry (Fig. [2]). In partic- 
ular, a transient cubic symmetric total density profile 
changes into an only rectangular symmetric total density 
profile. For a model of continuous orientations it can be 
expected to find an axial symmetry instead of the cubic 
symmetry. However, the spontaneous symmetry break- 
ing into a rectangular symmetric equilibrium state is ex- 
pected to be found for continuous orientational degrees 
of freedom, too. 

For an infinite system initially prepared within the 
isotropic-nematic two-phase coexistence region and one 
single laser beam of sufficient strength, instead of an 
array of laser beams, one would expect an equilibrium 
structure composed of three regions: A narrow high- 
density core at the position of the laser beam is sur- 
rounded by an annulus of nematic structure of total 
density (see Sec. IIIE[) which is surrounded by an 

annulus of isotropic structure of total density (see 
Sec. Ill E[) ; the sizes of the isotropic and nematic annuli 
are related to the initial state by means of the "lever 
rule" . In order to form this structure the system has to 
be large enough such that the amount of fiuid within the 
high-density core and the isotropic-nematic interface is 
negligible as compared to the total amount of fluid. How- 
ever, the two-dimensional grid applied within this work 
does not allow for sufficiently large single beam systems. 
In contrast, the smallness of the periodicity of the system 
considered in this work leads to a rather strong dilution 
of the fiuid in the space between the laser beams. 

In classical nucleation theory one considers the tem- 
poral evolution of fiuctuation induced "droplets" of the 
stable (nematic) phase surrounded by the metastable 
(isotropic) phase, where the densities are chosen as the 
two-phase coexistence values. The smallness of the peri- 
odicity of the current system, however, leads to a high- 
density core whose size equals the size of the laser beam 
and whose density is determined by the strength of the 
external potential. Therefore the high-density cores in 
this work are very different from the "droplets" consid- 
ered in classical nucleation theory. Moreover, the relax- 
ation paths here do not cross the " critical nucleus" , which 
corresponds to the minimal perturbation that lead to the 
equilibrium state and which is located at a saddle point 
of the free energy [53 |. Thus the current results can be 
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partly looked at in terms of the terminology of classical 
nucleation theory, because both are based on a free en- 
ergy function(al), however, the setup considered here is 
different from that of classical nucleation theory. 

After attaining a state close to the equilibrium state 
in the presence of the laser potential the laser beams are 
switched off and the relaxation of the (now) free fluid 
to equilibrium is traced. For sufhciently large initial su- 
persaturation of the initial state, one finds the central 
region of increased particle density stretching (Fig. [3]). 
The shape of these cores resemble factoids which have 
been found in dispersions of rodlike colloidal particles 
[53| . If the moving density fronts hit the periodic bound- 
aries of the modeled part of the system they merge with 
their "images" into a superstructure (Fig. [4|). As the 
nematic directors within the high-density cores of the in- 
dividual laser beams may point into different directions 
the merging of two neighboring cores is expected to give 
rise to one more relaxation process in which the merging 
cores reorient their colloidal particles along a common 
nematic director. This effect, however, is beyond the 



present study. 

In summary, a novel dynamic density functional the- 
ory for a model of fluids of platelike colloidal particles has 
been proposed and applied to the relaxation dynamics of 
the fluid under the influence of the potential of a two- 
dimensional array of laser beams. A rich phenomenology 
including the occurrence of individual approximately el- 
liptical high-density cores and the possible formation of 
complex superstructures has been found. 
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